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Abstract. Bell's theorem prevents local Kolmogorov-simulations of the singlet state 
of two spin-1/2 particles. We derive a positive lower bound for the L 2 -distance between 
the quantum mechanical spin singlet anticorrelation function cos and any of its classical 
approximants C formed by the stationary autocorrelation functions of mean-square- 
£h ' continuous, 27r-periodic, ±l-valued, stochastic processes. This bound is given by 

||C-cos|| > (1- 4f) /v^ps 0.133 95. 

PACS numbers: 03.65Ud 

> 

O ■ 1. Introduction 



Consider a Stern-Gerlach correlation experiment performed on the two distant members 
of a spin-singlet state of two spin 1/2 particles. The two apparatuses' magnetic 
fields point into the directions a, and b. Quantum mechanics tells that the stochastic 
correlation — Q (a, b) between the two measurements' results ±1 is given by —a ■ b, the 
negative scalar product between the two analyzing directions' unit vectors a and b. To 



get rid of the minus sign we employ in the following the anticorrelation function Q. 

Bell demonstrated that any local Kolmogorov simulation of such experiments is at 
variance with Q (a,b) = a ■ b. [1] The basic reason for this being his inequality, which 
subjects any triple of values C (a, b) , C (a, c) , and C (b, c) of the autocorrelation function 
C : S 2 x S 2 — > [—1, 1] , (a, b) i — y (fafb) of any local Kolmogorov simulation to 

\C(a,b) -C{a,c)\ <1-C{b,c). 

The quantum mechanical function Q, however, does not obey Bell's inequality for all 
choices of directions a, b, c. Thus Q = C cannot hold. 

This raises the question for how close the quantum mechanical function Q can 
be approximated by a local Kolmogorov autocorrelation function C. Clearly a notion 
of "closeness" is needed in order to make the question meaningful. The mean square 
deviation between C and Q averaged over a suitable set of pairs of directions is a natural 
choice, which we adopt. How should one choose the set of directions to be averaged over? 
Rotation invariance suggests to keep one point a G S 2 fixed and to average over all b 
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from a great circle or equivalently a half great circle through a, because from the SO (3) 
invariance3 of C it follows that C is completely determined by the mapping b t— y C (a, b) 
for any fixed a and all directions b lying on a meridian emanating from a. After all C 
is represented by a function of the (unoriented) angle between a and b only, i.e., there 
exists a function C : [0, it] — > [—1, 1] such that C (a, b) = C (9 a ,b) for all a, 5 G S 2 . Here 
6 a h G [0, 7r] is uniquely determined by a • b = cos a) &. 

Within the class of local Kolmogorov autocorrelation functions which are rotation 
invariant and continuous we have found that the mean square deviation of C from its 
quantum counterpart cos obeys the estimate 

J Q \C(8) -cos6\ 2 d6 > (l - /V2 « 0.133 95. 

The search for such "integrated" conditions which constrain the local Kolmogorov 
approximants C of the quantum anticorrelation function cos seems to have been initiated 
by Zukowski, who in [2] derived an upper bound for the modulus of the L 2 -scalar product 
between cos and C. See estimate (10) in [2], where it was taken for granted that certain 
integrability conditions are obeyed by the stochastic variables of the local Kolmogorov 
simulation under consideration. In [3] the idea has been generalized further in order to 
cover more general, entangled many particle states. 

In a first step we prove a scalar product bound, which in certain cases is equivalent 
to Zukowski's one. Our bound only holds for the singlet state but it does not rely on 
Zukowski's integrability assumption. The class of Kolmogorov simulations for which it 
holds is specified precisely in our lemma 5. From this we then derive in proposition 
6 a lower bound for the mean square distance between C and cos . Although fairly 
straightforward this latter bound seems to have remained unnoticed till now. It could 
serve the following purpose. 

Usually the decision upon the possibility of a local Kolmogorov simulation of a 
given set of spin singlett correlation measurements typically is based on the degree of 
violation of the Clausner-Horne-Shimony-Holt inequality for a finite set of analyzing 
field directions. [I] Our inequality might give a more sensitive criterion as to whether 
an experimentally determined singlet correlation function admits a local Kolmogorov 
simulation, since it involves the sum of (squared) deviations at a sufficient dense set of 
analyzing directions. If the mean square distance of the experimental data points from 
the quantum mechanically predicted values is estimated within a certain confidence level 
to be less than 0.13395, then a local Kolmogorov simulation can be ruled out within 
the same level of confidence. And this may well be the case even if all the individual 
vioaltions of the CHSH-inequality do not rule out locality with sufficient confidence. 

| This means that C (Ra, Rb) = C (a, b) holds for all directions a, b and for all R £ SO (3) . 
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2. Local Kolmogorov-simulation 

Let (x+iX-) denote an orthonormal basis of C 2 . Then the vector 

X = -j= (x+ ® x- - X- ® x+) 

represents the spin-singlet state of two spin-1/2 particles. A Stern-Gerlach experiment 
on each of the two constituents of the singlet state, which are assumed to be 
approximately localized in well separated regions, measures the two observables A = 
a (a) ® id and B = id® a (b) where a, b G § 2 C M 3 are two arbitrary vectors from the 
unit sphere S 2 . The unit vectors a, b specify the Stern-Gerlach apparatuses' orientations. 
For a = (a 1 , a 2 , a 3 ) G M 3 holds 



a (a) 



a 3 a 1 — ia 2 
a 1 + ia 2 —a 3 



The possible results of such measurements are the pairs (e, 77) G {1, —1} x {1, —1} . 
The quantum mechanically determined probability of outcome (e, rj) is given in terms 
of the euclidean scalar product a ■ b between a and b through 

l-erja-b 

Pa,b{e,v) = ^ • (1) 

This probability does not depend on the time order of the two Stern-Gerlach experiments 
since [A, B] = 0. 

Thus any choice of two directions a, b determines a probability measure on the 
event space {1,-1} x {1,-1}. The situation is analogous to a random experiment in 
which a pair of widely separated coins is tossed, while each coin of the pair is exposed 
to a magnetic field of direction a and b respectively. And the magnetic fields take an 
influence on the distribution of the possible results. 

The family of probability functions {p a ,b '■ Q>, b G S 2 } can be obtained as the 
distributions of stochastic variables {X a ^ : a, b G S 2 } 

X a ,b = {fa,b, 9a,b) : n ->• {1, "I} X {1, -1} 

on a single Kolmogorov probability space (O, W) as follows. [5] Take as the space of 
events Q the square [0, 1] x [0, 1] with the uniform distribution W. Then decompose the 
square Q into a set of four nonoverlapping rectangles, i.e. 

such that the area of the rectangle R a ^ (e, rj) has the value p a ,b (e, v) ■ Finally assume 

X a ,b (w) = (e, v) for w £ ^a,fe (e, v) ■ 

Then, by construction we obviously have 

Pa,b (e, 77) = W ({w G n : A ajfc (w) = (e, rj)}) . 

In such a model the outcome of each coin flip is determined by a randomly chosen 
"hidden variable" u G £1 in a way that in general the outcome also depends on both of the 
magnetic fields a and 6! How come that a coin gets influenced by faraway circumstances? 
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Trying to save locality, Bell considered the question whether there exists a 
probability space (O, W) with a stochastic variable 

X a , 6 =(/ a ,0 6 ):fi-).{l,-l}x{l,-l} 

for each pair (a, b) G § 2 x § 2 such that 

Pa,i, (e, 77) = W ({w G D, : / a (w) = e and # 6 (w) = 7?}) for all a, b G § 2 . (2) 

Such structure would constitute a /oca/ Kolraogorov simulation of the probability 
distributions generated by composite Stern-Gerlach experiments on a spin singlet state. 
The assumption that the stochastic variables f a ,gb depend on their respective Stern- 
Gerlach orientation a or 6 only is made in order to take care of the principle of locality. 

Now the equations ([I]) and (|2J) imply 

(fa9b) = 5^ £r]Pa ' b ( £ > V) = ~a-b 
e,ve{i,-i} 

for all a, b G S 2 . Yet Bell's theorem PQ rules out exactly that (f a 9b) — — a, • for all 
a, 6 G S 2 . Therefore a local Kolmogorov simulation of the singlet state does not exist. 
For the sake of completeness we spell out Bell's theorem precisely. 

Theorem 1 (Bell) Let (Q, W) be a probability space with two stochastic variables 
f a ,g a : Q {1,-1} for every a G S 2 . Then there exist points a, b G § 2 such that 
(fa9b) ^-a-b. 

Proof. That (f a gb) — —a-b cannot hold for all a, b G S 2 may be proven by contradiction. 
Choosing b = a the equation {f a g a ) = — 1 implies g a = —f a in the sense of stochastic 
variables, i.e., almost everywhere on Q. Thus, assuming (f a gb) — ~ a ■ b for all a, b G § 2 
leads to 

- (/oflb) = (fafb) =a-b 
for all a, 6 G S 2 . Bell noticed that, due to / 2 = 1, there holds 

|(/a/ 6 ) - (/a/c}| = |</a/ 6 (1 - / 6 /c)}| < 1 " (fbfc) 

for all a, 6, c G S 2 . Bell's famous inequality, however, is in contradiction with (f a fb) = a-b. 
Choose, e.g., three coplanar vectors a, 6, c with a • 6 = = b ■ c and a ■ c = 0. Then 

1 I A 1^1 a v^" 1 
= a ■ o — a • c < 1 — b ■ c = = — 

and thus the falsity 2 < v2 follows from Bell's inequality. ■ 
3. Quality of a classical singlet model 

Bell's theorem poses the following problem: Determine the infimum of the set of numbers 
77~T2 / I (fa9b) ~ [~a ■ b]\ 2 da do, 

(Alt) J§2 X§ 2 
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obtained from all local Kolmogorov models of a spin singlet state. Any such modej§| 
consists of a probability space (Q, W) and two families of {1, —1}- valued stochastic 
variables {f a : a G § 2 } and {g a : a G § 2 } such that (f a g a ) = — 1 for all a G S 2 . Thus 
again g a = — f a holds for all a G § 2 . Here da and denote the rotation invariant area 
element on the unit sphere normalized to Air. 

The solution to this problem would quantify and limit the optimal approximation 
to the quantum mechanical covariance function 

§ 2 x § 2 3 (a, b) ^ (x, a (a) <g> a (b) x ) = -a ■ b 

through classical singlet models. 

In this paper we address a somewhat simpler but related problem. We first confine 
the admissible direction vectors a, b from S 2 to a great circle S 1 C S 2 . Then we restrict 
to such ^-parametrized families of stochastic variables f a = —g a for all a G S 1 , for 
which there exists a continuous function C : R — > R such that 

(fafb) = C(a-b) holds for all a, b G S 1 . 

The assumption that (f a fb) depends on the scalar product a ■ b only amounts 
to postulating O (2)-invariance for the mapping (a, b) \-t (f a fb) , i-e., the relation 
(fiia fab) = (fafb) f° r eac h orthogonal mapping R which stabilizes the great circle S 1 . 
Under these premises we shall derive a positive lower bound for 



^T2 / Kfafb) - a ■ b\ 2 da db. 

27T) JgixS 1 



(2nY 

Here da and db denote the rotation invariant line element on the unit circle normalized 
to 2tt. 

An equivalent but simpler formulation of our problem is obtained by periodically 
parametrizing the circle S 1 through real numbers s and t, e.g., such that a (s) = 
(cos (s) , sin (s) , 0) and b (s) = (cos (t) , sin (t) , 0) . Then there exists a continuous 2n- 
periodic function C : M — > R with C (a ■ b) = C (t — s) and we obtain 



(27T) 2 



/ Kfafb) ~a-b\ 2 da db 



(27T) 2 



2tt /.2tt 

Jo 



C(t-s) - cos (t-s)| 2 ds dt 



r lew- cos (t)\ 2 dt. 

2tt Jo 



|C(t) - cos (t)| 2 dt, 



Thus we try to find a positive lower bound for 
1 Z" 2 ^ 
2^7o 

where the continuous function C : R — )• R is related to a 27r-periodic, {1, — l}-valued, 
stochastic process {f s : s G R} through C (i) = (f s f s+t ) for all s, t G R. 



§ We denote such models as classical singlet models. 
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Notice that, because of 

C(-t) = (f s f s - t ) = (f S -J S - t+t ) = (f s f s+t ) =C(t), 

the function C is even and it also obeys C (0) = 1. 

If, in addition to C (t) = (f s f s+t ) , the mapping s (f s ) is constant, the process 
{f s : s £ 1} is called stationary in the wide sense. [6], [7] In case of (f s ) = for all s G R, 
the function C specializes to the autocorrelation function of a wide-sense-stationary 
process {f s : s G R} . However we will not need any assumption on the expectation 
values (f s ) . 

The following results from the theory of stationary processes make it clear that 
insistence on the continuity of C is much less a restriction than it appears to be. They 
show that the condition of continuity of C can be replaced by the seemingly weaker 
condition that C is continuous at 0. See, e.g., sect. 8.10 from [6]. For the sake of 
completeness we include the proofs. 

Lemma 2 Let {f s : s G R} be a {1, — 1} -valued stochastic process such that there exists 
a function C : R — > R with C (t) = (f s f s +t) f or a ll s, t G R. Then C is continuous 
everywhere if and only if it is continuous at 0. 

Proof. Observe first that 

C(t + e)-C(t) = (f f t+E ) - (f f t ) = (f (f t+£ - f t )) . 
Thus we have, due to the Cauchy-Schwarz inequality and due to ft = 1, that 

\C(t + e)-C (t)\ 2 = \(f (f t+£ - f t ))\ 2 < </ 2 > ((f t+£ - f t f) 

— ({ft+e — ft) ) = (ft+e + ft ~ %ftft+s) 

= 2(1- C7(e)) = 2(C7(0)- (7(e)). 

From this it follows that 

lim \C{t + e)-C (t)\ < /21im (C(0)-C( £ )). 

Thus \im £ ^ \C (t + e) — C (t)\ = if lim e _> C (s) = C (0) . Clearly, if C is continuous 
everywhere it is in particular continuous at 0. ■ 

Definition 3 A stochastic process {f s : s G R} , with lim e _> ((ft+e — ft) 2 ) = for all 
t G R is called mean- square- continuous. 

Lemma 4 A stochastic process {f s : s G R} with C (t) = (f s f s+t ) for all s, t G R with 
values in {1,-1} is mean-square-continuous if and only if C is continuous at 0. 

Proof. In the case of a {1, — l}-valued process with stationary correlation function C 
holds 

lim<(/, +£ -/,) 2 > =Mj(/^ + {/ i V2(/ H J t )]=21im(l-C(e)) 
= 21im(C(0)-C(e)). 
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Thus in the case of a {1, —1}- valued stochastic process {/ s : s G R} , for which 
there exists a function C : R — > R such that C (t) = (f s f s+t ) holds for all s,t G R, the 
following three conditions are equivalent: 

(i) C is continuous. 

(ii) C is continuous at 0. 

(iii) {f s : s G R} is mean-square-continuous. 



4. Lower bound for the quality 

Our main tool is the following estimate for the Fourier coefficients of the correlation 
functions C which appear in the present context. 

Lemma 5 Let {f s : s G R} be a 2n -periodic, {1, —1} -valued, mean-square- continuous 
stochastic process such that C (t) = (f s f s+t ) holds for all s, t G R. Then all the Fourier 
coefficients Ck of C exist and for all k G Z the following estimates hold 



2tt 



— / cos (kt) C (t) dt < (-} . 



0<c k = — e ' M C(t)dt = — 
2tt Jo 

Proof. Since the function C is continuous, its Fourier coefficients exist. Since C is real 
valued and even, there holds c& = c^k = c&. Thus the mapping k >->■ Ck is real valued 
and even too. In particular because of Ck G R we have 

r 2TT 



f 

Jo 



e- lkt C (t) dt 



e-' lkt C (t) dt = 
cos (kt) C (t) dt 



(cos (kt) - i sin (kt)) C (t) dt 



2tt 



cos (kt) C (t) dt. 



Now for the estimate c k > : 



2vr/ 
1 

2^ 



e- M (f s f s+t )dt 



2ir 



2tt 



-ikt 



1 



2tt 



1 

2^ 



27T 



2tt 



iks 



\ks 



2tt 



2tt Jo 

(/ s / s+t ) d S 



(/,/.+*> ds dt 



o 

N N 



2tt 



e~ lkr (f s fr) dr ds 



m=l ri=l 

TV 



lim 



2^ ^ 



ri=l 



iV 



> 0. 



The last line's sum can be bounded from above as follows. First observe that for 
all r G R it holds 



Eexp (-ik2n%) 
hn- 



71=1 



N 



kr ™ exp (-ifcgTrf) 



n=l 



N 
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Now there exists a number r N e [0, ^) depending on N, such that 

k " exp (-ifegrg) 
2^ - 72^ 



n=l 



iV 



n=l 



iV 



holds. It then follows for such r/y that 

Eexp (-ifc27rf ) 
^ i2^i 



n=l 



AT 



N 

E 

n=l 



cos(fc(27rf -Tjy)) 



Since /W^- assumes values from {—1, 1} only, we obtain from the triangle inequality 



N 

E 



cos (k (27rf - rjy)) 



1 



N 



2tt|cos (fc (27rf -rjv)) 
iV 



n=l n=l 

The points Z]y = {2ttj^ — tn ■ n = 1, . . . iV"} partition the interval [0, 2tx]—tn whose 
length is 2-7T. The mesh 2n/N of Zn tends to zero for iV — > oo. Furthermore 2tt is a 
periode of the function x i— )■ |cos . Therefore the sum on the right hand side of this 
last inequality converges towards the Riemannian integral 

1 f 2lT 

— / \cos(kt)\dt 
27r Jo 

for A" — > oo. For each e > there thus exists a number A/q, such that for all N > Nq 



N 

E 



cos 



71=1 

For G Z \ we have 



(k2n%-T N )\ < 1 
A 7 " ~ 2tt 



2tt 



| cos (/ct)| dt 



cos (Art) | dt = 4A: / cos (fct) dt = 4 sin (kt) \^ /2k = 
o Jo 

Thus for any e > there exists a No, such that for all N > Nq holds 



V"* exp (- ik2n Jr ) * 

^ hir% 



n=l 



N 



<-+s. 

7T 



We therefore have proven for any fceZ\0 that 



lim 

N^oo 



ex P (- i/c27r f ) f 

2^ ^ 



n=l 



AT 



< 



For k = the estimate Co < (-) follows from 
1 /" 27r 

■ 

Note that in this proof we did not interchange the limiting process of integration 
with the probabilistic expectation value, which in general also involves a limit process. 
Such an interchange can be misleading since the realizations t (->■ f t (ou) need not be 
integrable for almost all u> G Q. If, however, the two limits can be interchanged, the 
proof gets abbreviated considerably. [2] , [B] 
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From the lemma's estimate for the case k = 1 namely < C\ < (^) 2 we now 
obtain our lower bound for the L 2 -distance between the quantum mechanical spin singlet 
correlation function and its classical approximants. 

Proposition 6 Let {f s : s G K} be a 2n -periodic, {1, —1} -valued, mean- square- con- 
tinuous stochastic process such that C (t) = (f s f s+ t) holds for all s,t e R. Then the 
mean square deviation of C from the quantum mechanical correlation function cos obeys 



lie- 

Proof. Note that 



cos 



Lj\ cit) . coe{t) f At> _\zA 



i r 2w i 

= -J co S *( t )d t = 5 . 



0.133 95. 



cos 



We thus can decompose C into a component parallel to cos and one orthogonal to it 
according to 

C = V2 cos (y/2 cos, &j + (c - v^cos (y/2 cos, C 

Here the scalar product between two continuous functions f,g:] 
of 2-7T is denoted by 



— >■ C with a periode 



(f,g) = ^JJ f(t)g(t)dt. 



From the estimate 



Cl 



1 

2^ 



27T 4 

cos (t) C (t) dt< — 

TV 2 



it follows that 



\\C - cosf = — / (C (*) - cos (t)) 2 dt 
2tt Jo 

cos, C ) 
cos ( \Pl cos, C ) — COS 



> 



V^cos/v 7 ^ cos, C) — COS 



;i-2(cos,C)) 2 ||cos| 



= i(l-2(cos,C)) 2 . 
Because of < 2 (cos, C) = 2ci < ^ = 0.810 57 we finally have 



|C — cos 1 1 > 



1 - -2- 

v/2 



0.133 95. 



From this proof it is obvious that the estimate \\C — cos|| > (l — %) /y/2 is 
saturated if and only if C is proportional to cos which, because of C (0) = 1, in turn 
implies C = cos . Thus the estimate cannot be saturated and stronger estimates might 
exist. 
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5. Bell's example 



Bell [1] constructed a local spin singlet model with the 27r-periodic autocorrelation 
function given by C : R — > [—1, 1]: 

C(t) = 1 -2^1 for < Itl < 7i. 

C is continuous and even. Bell's stochastic variables {f a : a G § 2 } are defined on the set 
S 2 endowed with the uniform distribution. They are given through 



/«(") 



1 for a; • a > 
-1 otherwise 



and indeed yield (f a fb) = 1 — 2^ with 9 G [0, tt] such that a ■ b = cosO holds. 
The Fourier coefficient of C is given for k G Z \ through 



k 2tt 



-ifct 



7T 



2tt 



1 - 2^ dt = — / cos (Art) 1 - 2-^ )dt 



7T 



- r cos (Art) f 1 - 2- J dt = r 4- [s in (**)] fl - 2- J dt 



1 

Ar7r 



sin (Art) 1-2 



7T 



+ 



7T 



sin (Art) dt 



2 2 
—2 cos (Art) |o = -—-2 



1 -(-!)' 



7-^7 for odd A; 

(KTC) 

for even /c 

Obviously Co = holds. The mapping k (-)■ indeed is real valued and even. The 
estimate < Ck < 4/n 2 is realized and in the case of ci saturated. Therefore it holds 
that 

\\C - cosf = \\Cf + ||cosf - 2 (cos, C) 
,2 , 1 



2 
1 

2 



2c 



|C|| 2 + ~ - -j. 

7T 



The value of ||C|| is given by 

,2 1 



c 



2tt 
1 

7T 
1 



C 2 (t) dt 



7T 



C 2 (t) dt 



1-2-1 dt 

7T 



2ar) 2 dx 



Thus we have 



y 2 dy 



1 1 



y 2 dy = 77- 
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and in consequence 

IIC-cosI 



- = 0.150 88. 

6 7T 



Note that C has the following particularly simple uniformly converging Fourier 
series representation 



C (t) = ( c ^ kt + ^- k e- m ) = c * + e~ ikt ) 

k=l k=l 

oo oo 

= 2cfc cos (kt) = 2c 2/ t + i cos ((2k + 1) t) 



k=l k=0 

o OO 



7T* 



E 

fc=0 



cos ((2Jfe + 1) t) 



(2k + 1) 



cos (t) + 



cos (3t) cos (5t) 



+ 



+ . . . 
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